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ABSTF  ACT.  The  value  of  large  software  products  can  be  cheaply 
increased  by  adding  restyled  interfaces  that  attract  new  users. 
As  axaaples  of  this  approach,  a set  of  graphics  primitives  and  a 
language  precompiler  for  scientific  computation  are  described. 
These  two  systems  include  a general  user-defined  coordinate 
system  instead  of  numerous  system  settings,  indention  to  specify 
block  structure,  a modified  indexing  convention  for  array 
panmeters,  a syntax  for  n-and-a-half-times-' round  loops,  and 
engineering  format  for  real  constants;  most  of  all,  they  strive 
to  be  as  small  as  possible. 


0.3  PHILOSOPHY.  Kernighan  and  Plauger  f 1976  ] describe 
explicitly  and  by  example  three  precepts  cf  the  Software  Tools 
philosophy: 

trim  out  the  inessentials 
build  it  adaptively 
- let  someone  else  do  the  hard  part 
Two  more  examples,  driven  by  the  same  philosophy,  are  given 
below.  The  basic  idea  is  to  obtain  high  leverage  by  taking  an 
existing,  powerful  piece  of  software  and  make  it  useful  to  more 
people  by  designing  a new  interface.  Webster's  calls  this 
process  facelifting:  "a  restyling  intended  to  increase  comfort 
or  salability." 


1.0  JUSTIFICATION  FOP  STILL  ANOTHER  PROGRAMMING  LANGUAGE. 
Fortran  will  no  doubt  remain  for  many  years  the  most  important 
programming  language  for  scientific  computation.  Nhen  used 
carefully  and  with  discipline,  it  yields  remarkably  portable 
codes;  this  is  its  greatest  virtue.  But,  as  programmers  have 
coaplained  for  years,  it  also  has  many  faults: 

- awkward  syntax  for  stateaeats,  strings,  names 

- primitive  control  structures 

- DO  loop  restrictions 

- no  macros 

Fortran  preprocessors,  such  as  MORTRAN  ( Cock*  Shustek  1975],  have 
eliainated  many  of  these  disadvantages  and  therefore  have  become 
very  popular.  Unfortunately,  they  reduce  portability  soaewhat, 
since  either  the  preprocessor  oust  be  Installed  at  the  new  site 


o 


or  illegible  'object'  Fortran  sent  there.  Hore  iaportantly,  such 
preprocessor s have  only  a sinor  effect  on  inherent  probleas  of 
Fortran: 

- dynaaic  allocation  is  either  unavailable  or  requires  the  use 
of  rather  confusing  tricks 

- no  PH7CID0PE  VARIABLE  type 

- no  STHUCTUFE  type 

(Labelled  cosaon  blocks,  since  they  do  not  use  the 
ccabinatorial  possibilities  of  procedure  paras  et  er  izat  ion, 
are  less  flexible.) 

- no  0-origin  indexing 

- array  bound  infor nation  is  not  aotoaatically  passed 

- no  vector  operations 

- no  recursion 

The  PORT  library  Bakes  dynaaic  allocation  one  of  its  aost 
advertised  features:  "He  have  found  that  use  of  dynaaic  storage 
allocation  in  POBT  leads  to  aore  clearly  structured  prograas, 
cleaner  calling  sequences,  iaproved  aeaory  utilization,  and 
better  error  detection."  [ Fox*Hall+Schryer  1977]  Adding  a stack 
tc  Portran  is  a nessy  affair,  however,  as  shown  in  figure  1, 
which  contains  two  alternate  nethods  in  PCFT  for  allocating  an 


SUBROUTINE  LBB(A,N) 


COMMON  /CSTAK/DSTAK  ( 500) 

DOUBLE  PRECISION  DSTAK 
INTEGER  ISTAK(IOOO) 

REAL  A(l) 

REAL  RSTAKU000) 

EQUIVALENCE  (DSTAK(I),ISTAK(D) 
EQUIVALENCE  (DSTAK(I).RSTAKO)) 

II  = ISTKGT(2*N,2) 

IR  = ISTKGT(N,3) 


SUBROUTINE  l.BB(A,N) 

COMMON  /CST A K / DST A K ( 500) 

DOUBLE  PRECISION  DSTAK 
INTEGER  ISTACKOOOO) 

REAL  A(l) 

REAL  RSTAKOOOO) 

EQUIVALENCE  (DSTAKG).ISTAK(I)) 
EQUIVALENCE  (DSTAK (l),RSTAK(l)) 


( code  referring  to  RSTAKdR  + n)  and  ISTAKdl  + m) 
probably  ending  with  code  to  store  the  stuff 
from  the  real  scratch  storage  into  array  A ) 


II  - ISTKGT(2*N,2) 

IR  = ISTKGT(N,3) 

CALL  LIBB(A,ISTAKdl),RSTAK(IR),N) 


CALL  ISTKRL(2) 

RETURN 

END 


CALL  ISTKRL.(2) 

RETURN 

END 


figure  1 


2 


INTEGER  and  BEAL  array. 

Other  proposals  are  even  aore  coaplicated.  (After  a 7 page 
description  of  DYNOSCP,  Huybrech ts(  1977  ] states:  "This  paper  < 

qives  cnly  the  basic  features  of  the  PYNOSOB  systea.  A aore 
sophisticated  use  allows  the  user,  once  he  Is  faailiarized  with 
the  systea,  to  iaprove  greatly  the  speed  of  prograas  using  it.") 

PL/I,  which  is  now  becoaing  fairly  widely  available  in  soae  fora, 
overcoaes  all  these  difficulties.  However,  so  huge  a language 
tends  to  overwhela  people,  and  because  of  tricky  precision  rules, 
silant  type  conversions  (as  in  I*J=0;l,  and  the  like,  learning 
only  part  of  the  language  is  dangerous. 

other  languages,  while  beautifully  designed,  have  their  own 
flaws.  For  exaaple,  Algol  H does  not  have  a robust  interface  to 
Pcrtran;  in  addition  to  this  (Fohilner  1977  ],  Pascal  places 
painful  restrictions  on  arrays.  < 


1.1  T.  Thus  another  apprcach  seeas  warranted,  which  can  coabine 
the  needed  features  of  PL/I,  the  deliberate  syntax  of  ALGOL,  and 
the  low  inpleaentation  cost  of  the  Fortran  preprocessors.  Such 
an  approach  has  produced  the  language  T,  intended  to  assist  in 
the  inpleaentation  and  docu  sent  a tion  of  algorithas  for  scientific 
coaputation.  The  principal  aias  have  been  ease  of  reading  and 
writing,  low  inpleaentation  cost,  and  reasonable  efficiency. 

Appendix?  gives  the  foraal  language  proposal,  specifying  the 
syntax  according  to  Mirth’s  proposal  [1977  ],  Since  T is  siailar 
to  Fortran,  Algol  60,  and  PL/I,  a coaplete  specification  of  the 
seadntics  nay  be  oaltted  without  confusion.  To  provide  the 
heuristics  behind  the  design  choices  and  to  give  an  overview  of 
the  language,  various  aspects  of  the  folowing  exaaple  will  be 
discussed. 


TRIP  EAR 

• exaaple  of  T and  G systeas; 

• various  views  of  the  sua  of  three  Gaussian  peaks; 
t Eric  Grosse  Stanford  University 

REAL:  AZIN,  ELEV,  • VIEMI RG  ANGLES  FOR  SURFACE  PLOT 
BELEFB , A ESEFR  , t ERROR  TOLERANCES  FOR  ODE 

• IRCBPERCZNT  VARIABLES  OF  TRAJECTORY 

• 2 NORN  OF  THE  GRADIENT 

• CORNERS  OF  RECTANGULAR  DOHA  IN  OF  FUNCTION 

• FOCAL  POINT  FCR  SURFACE  PLOT 

• COORDINATE  TRANSFORMATION  PARAMETERS 

• LOCATION  AND  GRADIENT  FOR  TRAJECTORY 

I NTEGE  F ( 5)  : ODEIWOFR 


T,  TOUT, 

NOB  HYP 

REAL  (2)  : LL,  UR  , 
ORIGIN  , 

XO,  SCALE, 
Y,  YP 

RE  AL  ( 142)  : ODEHOPK 
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OEFINE  (P, 20)  • density  of  P samples; 

REAL(-P:  P,-P:  P)  : P TABLE 

9 EAL  ( 3)  : LEVEL  • CONTOOB  LEVELS  > 

INTEGER:  I,  J, 

IPLAG  * DIAGNOSTICS  FLAG  FOB  ODE 

STRUCTURE:  PARAH  • LOCATIONS,  HEIGHTS,  AND  WIDTHS  OF  PEAKS 

REAL (3,2):  X 
REAL  (3):  H,  W 

STRUCTURE:  PF  • PLOT  FILE 

INTEGE  B (500)  : WORK 

PROCEDURE:  GOPEN,  GCLOSE,  GPICT,  GCONT , GSOBF,  GLTTPE, 

GJUNP,  GDRAH,  GTRAN1 
PORTRAN  PROCEDURE:  CDE,  DP,  STASH 
PROCEDURE  ()  REAL:  F 

# SET  UP  PARAMETERS 
BLANK  SEPARATION  (2) 

REAL  DIGITS  (3) 

GET  DATA  (A ZI H ,E  LEV ) 

PUT  DATA  ( AZ  IF.  , ELEV) 

X (1,1)  :*  o 

X (1,2)  :*  0.5 
X ( 2 , 1)  :=  -o.  4 3 30  1 2 7 0 2 
X (2,2)  :*  -0.25 
X ( 3,  1)  :=  “T  (2,1) 

X (3,2)  :=  X (2,  2) 

PUT  DATA  ARRAY  (X) 

GET  ARRAY  (H) 

PUT  DATA  ARRAY  (H) 

GET  ARRAY  (H) 

PUT  DATA  ARBAY  (W ) 

STASH  (X,H,H) 

FOR  ( -P  <=  I <-  P ) 

Y (1)  : = FLOAT  (I)  / P 
FOR  ( - P <*  J <*  P ) 

Y (2)  :=  FLOAT  (J)  / P 
F T ABLE  (I , J)  :*  F ( Y,  PARA  H) 


• SURFACE  PLOT 
GOPEN  (*  VEP12FF*  ,PF) 

3 PICT  ( PF) 

LL  :*  -1 
UR  :=  1 
ORIGIN  :*  0.5 

GSURF  (LL,  U R,  FT  ABLE,  AZIH  , ELEV  , ORIGIN,  0. 25,  PF) 


T- 


• CONTOUR  PLOT 
GPICT(PP) 

SCALE  :*  0.3333 

X0  :*  -0.  5/SCALE  (1) 

GTBAN1  (XO,  SCALE,  PP) 

GET  AF8AT  (LEVEL) 

PUT  DATA  ABBAT (LEVEL) 

SCOHT  (LL,UR  ,PTAB LE, LEVEL,  PP) 

GLTTPE  ( ' DOT*  ,PP) 

GET  AFBAT  (LEVEL) 

PUT  DATA  ABBAT  (LEVEL) 

SCOUT  (LL.tlB  .PTABLE,  LEVEL,  PP) 

• COMPUTE  AND  PLOT  TBA JICTORT 
BE1BBB  :*  10 (-6) 

GLTTPE  ('SOLID  ' , PP) 

ABSEBB  :*  10  (-6) 

WHILE ( - END  OP  INPUT  ) 

GET  ABBAT  ( T ) 

PUT  DATA  ABP  A T ( T ) 

T :«  0 
GJU  HP  (T,  PP) 

IPLAG  :«  1 

WHIL E ( NOBP  TP  > 1 (-  3)  6 1 <*  IPLAG  6 IPLAG<»3  ) 

TOUT  :«  T ♦ 10  (-3)  /NOBHTP 

CDE  (DP,  2,  T,  T,  TOUT, BEL  EBB  ,ABSERB  , IP  LAG  ,ODEHOBK,  ODEIHOBK) 
CASE 

2 « IPLAG 

GDB AN  (T,  PP) 

3 * IPLAG 

PUT  ('ODE  DECIDED  EIBOP  TOLERANCES  BEBE  TOO  SHALL.') 
PUT (' NEN  VALUES:') 

PUT  DATA  (BEL EBB, ABSEBB) 

ELSE 

PUT('OCE  RETURNED  THE  EBBCB  PLAG:') 

PUT  DATA  (IPLAG) 

PIBST 

DP(T,  T, TP) 

NOBHTP  :*  NOP H2 ( TP) 

GC  LOSE(PP) 

P ( T,  PABAH  ) Z 
BEAL  ()  : T 
REAL:  Z,  N0BH5Q 
STRUCTURE:  PABAH 
REAL  (3,2):  X 
REAL  (3)  : H,  W 
INTEGER:  I 
Z :*  0 

POB  ( 1 <■  I <*  3 ) 

NOBHSQ  :«  (T  (1)  -X(I,  1)  )**2  ♦ (T  (2) -X  (I,  2) ) **2 
Z :«  Z ♦ H (I ) ♦EXP  (-0 . 5*H ( I)  *NOR  HSQ) 
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1.2  CONTBJL  AND  CTHEP  SYNTAX.  Perhaps  the  lost  striking  feature 
the  Algol  veteran  sees  in  this  exaaple  is  the  ooaplete  absence  of 
BEGINS  and  ENDs.  Not  only  is  the  text  indented,  bat  the 
indention  actually  specifies  the  tlock  structure  of  the  prograa. 
Such  a scheae  was  apparently  first  proposed  by  Landin  ( 1966]. 
Except  for  an  endorseaent  by  Rnuth  [1974],  the  idea  seeas  to  have 
been  largely  ignored.. 

Ideally,  the  text  editor  would  recognize  tree-structured  prograas 
[Hansen  197  1].  In  practice,  text  editors  tend  to  be  line 
oriented  so  that  loving  lines  about  in  an  Indented  prograa 
requires  cuabersoae  aanipulation  of  leading  blanks.  Therefore 
the  current  iapleaentat  ion  of  T urns  BEG  T N and  END  lines, 
translating  to  indention  on  output.  Thus  the  input 
STB  fJCT  OPE:  PAHAE 
(( 

BEAL  (3 ,2)  : X 
BEAL  ( 3) : H,  tf 
) ) 

praduces  the  output 

STFU-TtJBE:  PABAP 
REAL  (3,2):  X 
PEAL  (3)  : H,  H 

whatever  the  iapleaentation , the  key  idea  is  to  force  the  block 
structure  and  the  indention  to  be  autoaatlcally  the  saae,  and  to 
reduce  clutter  froa  redundant  keywords. 

Blanks  are  Insignificant  outside  of  strings,  natheaatical  tables 
have  long  used  blanks  inside  nuaeric  constants,  as  in 
PI  :=  3.  14159  26535  89793 

for  readability.  Blanks  in  identifiers  alsc  can  iaprove 
readability,  while  reducing  the  chance  of  aisspelling  and  easing 
the  pain  of  naae  length  restrictions  iaposed  by  the  local 
operating  systea. 

In  accordance  with  the  recoaaen dations  of  ScowentWlchaann  [ 1973], 
coaaents  start  with  a special  character,  #,  and  run  to  the  end  of 
the  physical  line. 

The  saall  reserved  word  list  eliainates  the  need  for  a stropping 
convention.  The  psychological  advantages  of  this  approach  nave 
been  elaborated  by  Hansen  [ 1973  ]. 

The  fora  of  the  assianaent.  and  procedure  call  stateaents  follows 
the  clean,  clear  style  of  Algol  6C.  To  sake  macros  aore 
understandable,  their  syntax  and  seaantics  aatch  those  of 
procedures  as  closely  as  pcssible. 

In  addition  to  noraal  stateaent  sequencing  and  procedure  calls, 
three  control  structures  are  provided.  Tne  CASE  and  WHILE 
stateaents  are  illustrated  in  this  typical  prograa  aegaent: 


ia<i  tAnTi 


9 HZ  IE  ( BOB  HY  P > 1 (-3)  6 1<*IPLAG  t IPUG<*3  ) 

TOOT  : * T ♦ 10  (-3)/HOBHYP 

ODE  (OP,  2,  Y, T, TOUT ,RELEBF, IBS ERB,  IPL  AG,  ODE  WORK  , ODE  1 90B  K) 

CASE 

2 * IPLAG 

GDBAH  (Y,PP) 

3 * IPLAG 

POTCODE  DECIDED  ERBOB  TOIEBIBCES  9E BE  TOO  SHALL.') 

PUT  ( ' ME9  VALUES:') 

PUT  DATA  (BELERB,  ABSERB) 

ELSE 

PUT  ( 'CDE  FETOBH  ID  THE  EBBOB  P1AG:') 

PUT  DATA  (IPLAG  ) 

PIP  ST 
DP(T,Y,  YP) 

BOBPYP  :=  N0BH2  (TP) 

The  CASE  statement  is  modelled  after  the  conditional  expression 
of  LISP;  the  boolean  expressions  are  evaluated  in  sequence  until 
one  evaluates  to  YES,  or  until  ELSE  is  encountered.  The  use  of 
indention  aakes  it  easy  to  visually  find  the  relevant  boolean 
expression  and  the  end  of  the  statement. 

One  unusual  feature  of  the  WHILE  loops  is  the  optional  PIBST 
■ truer,  which  specifies  where  the  loop  is  to  be  entered.  In  the 
exaaple  above,  the  nora  of  the  gradient,  NOBHYP,  is  computed 
before  the  loop  test  is  evaluated.  Thus  the  loop  condition, 
which  often  provides  a valuable  hint  about  the  loop  invariant, 
appears  proainently  at  the  top  of  the  loop,  and  yet  the  coaaon  n- 
and-a -half-tiaes-' round  loop  can  still  be  easily  expressed. 

The  PCS  stateaent  adheres  as  closely  as  practical  to  coaaon 
aatheaatical  practice. 

PDF ( 1 <*  I <*  1 ) 

NOFHSQ  :*  (Y  (1)-X(I,  1)  )**2  ♦ (Y  (2) -X  (T,  2)  ) **2 
Z I*  Z ♦ H (I)  *E7P(-0. 5*9(1)  ♦HOBHSQ) 

Sevaral  years  experier.ee  with  these  control  constructs  has 
deaonstrated  then  to  be  adequately  efficient  and  such  easier  to 
naintain  than  the  alternatives. 

Procedure  nesting  is  not  used  for  two  reasons.  first,  textual 
nesting  that  extends  over  aany  pages  is  difficult  for  a huaan  to 
keep  track  of.  Second,  progress  typically  contain  several  high 
level  procedures  calling  a single  priaitivc,  so  a tree 
representation  is  inappropriate  anyway. 

By  reaoving  the  nesting  of  procedures,  however,  we  worsen  the 
problea  of  entry  point  hidinq  that  arises  when  coabinlng  progress 
froa  aany  sources  into  a sinqle  library.  A solution  to  this 
problea  is  to  have  an  official  Dane  for  each  procedure,  coded 
along  the  lines  of  IHSL,  and  also  a aore  aneaonic  nick  naae 
(which  users  can  pick  for  theaselves  if  they  like).  The  aacro 
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processor  which  is  built  into  T can  then  be  used  to  change  all 
occurences  of  the  nick  na  aes  into  the  corresponding  official 
na  ess. 


1.3  DSC  LA  RATION  S . The  fundaaental  scalar  types  are  INTEGER, 
REAL,  and  COMPLEX,  fro ■ which  arrays  and  structures  aay  be  bnilt 
up.  As  the  exaaple^ 

BEAL  (-P:P,-P:P) 

illustrates,  general  upper  and  lower  bounds  are  allowed. 

The  upper  bound  expression  is  oaitted  for  a foraal  array 
parameter,  so  that  an  appropriate  value  car  be  taken  froa  the 
length  of  the  corresponding  actual  array  arguaent.  The  origin  of 
an  actual  array  arguaent  need  not  aatch  the  origin  of  the 
ccrre  sponding  foraal  array  paraaeter.  For  example,  if  the  actual 
arguaent  A was  declared  FEAL(0:7):  A and  the  foraal  paraaeter  B 
was  declared  HE  AL  ()  : P , then  B(8)  will  correspond  to  A (7).  Host 
languages,  when  they  allow  lower  bounds  at  all,  do  not  perait 
this  flexibility,  which  is  used  in  the  exaaple  prograa  when  a 
aatrix  with  lower  bound  -P  is  passed  to  a general  purpose  library 
routine  which  assuaes  a lower  bound  of  0. 

Structures  of  arbitrary  depth  nay  be  declared.  As  the  exaaples 
STRUCTURE:  P ARAB 
BEA  L ( 3,  2)  : X 
B BA  L (3):  H,  W 
STRUCTURE:  PF 

INTEGER  (500)  : BOR  K 

suggest,  structures  are  useful  passing  collections  of  related 
data,  without  the  need  for  long  paraaeter  lists.  This  aakes 
feasible  the  prohibition  of  global  variables  in  a drastic  atteapt 
to  narrow  and  sake  aore  explicit  the  interface  between 
procedures.  Euclid  ( Popekvothers  1977  ] has  emphasized  the 
importance  of  visibility  of  names. 

The  graphics  procedures  which  use  the  BCBN  vector  of  the  exaaple 
are  able  to  divide  up  the  space  into  convenient  units.  This 
capability,  which  would  be  possible  in  PL/I  only  through  the  use 
of  pointers,  encourages  inforaation  hiding  and  abstraction. 

PPOCEDUBE  VARIABLES  allow  the  nanes  of  procedures  to  be  saved,  an 
essential  feature  for  applications  like  the  user-specified 
coordinate  transform ation  described  in  the  graphics  systea  below. 

The  importance  of  existing  Fortran  software  is  recognized  by 
providing  for  FORTRAN  PROCEDURES  as  an  integral  part  of  the 
language.  The  current  iapleeentation  of  T performs  this  linkage 
in  a aore  efficient  way  than  the  naive  user  of  PL/I  would  be 
likaly  to  discover. 
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A novel  syntax  is  introduced  for  function  returns.  Since 
procedures  say  be  recursive,  Fortran's  convention  of  using  the 
function  naie  as  variable  cannot  be  followed.  Instead,  the 
procedure  header  declares  a return  variable  lust  like  any  other 
paraseter: 

F ( Y,  PA  BA  H ) Z 
PEAK):  Y 
HEAL:  Z 


1.4  IHPDT/D OTPUT.  Beginners  often  find  Fortran's  input/output 
the  sost  difficult  part  of  the  language,  and  even  seasoned 
prograssers  are  teapted  to  lust  print  unlabelled  nusbers,  often 
tc  sore  digits  than  Justified  by  the  problea,  because  foraatting 
is  so  tedious.  PL/ 1' s list  and  data  directed  I/O  is  so  auch 
easier  to  use  that  it  was  wholeheartedly  adopted  in  T.  By 
providing  procedures  for  acdifying  the  nuaber  of  deciaal  places 
and  the  nuaber  of  separating  blanks  to  be  output,  no  edit-drected 
I/O  is  needed.  Special  stateaents  are  provided  for  array  I/O  so 
that,  unlike  FI./I*  arrays  can  be  printed  in  orderly  fashion 
without  explicit  foraatting. 

Since  alaost  as  auch  tiae  is  spent  in  scientific  coaputation 
staring  at  pages  cf  nuabers  as  at  pages  of  prograa  text,  auch 
thought  was  given  to  the  best  foraat  for  displaying  nuabers. 

Tn  jccordance  with  the  "engineering  foraat"  used  on  Hewlett- 
Packard  calculators  and  with  standard  aetric  practice  r GH  Service 
Section  1977],  exponents  are  forced  to  be  aultiples  of  3.  As 
figure  2,  an  excerpt  froa  the  exa  aple  prograa's  output,  shows, 
thi3  convention  has  a histograaa  ing  effect  that  concentrates  the 
inforaation  in  the  leading  digit,  as  opposed  to  splitting  it 
between  the  leading  digit  and  the  exponent,  which  are  often 
separated  by  14  coluaas.  The  use  of  parentheses  to  surround  the 
exponent,  like  the  legality  of  ia  tedded  blanks,  was  suggested  by 
aatheaatical  tables.  This  notation  separates  the  exponent  froa 
the  aantissa  acre  distinctly  than  the  usual  E foraat. 


1.5  DISCUSSION . 

Following  Kernighan*Plauger  [ 1976  ],  the  initial  lap  Indentation  is 
unsophisticated  [Cower  1978].  Nevertheless,  the  preprocessing  is 
les3  costly  than  the  FL/I  coapile,  so  the  cverall  results  are 
guite  satisfactory.  (The  evaluation  looks  even  better  if  one 
• ccapaces  PL/I  ♦ T against  PL/I  ♦ PL/I's  aacro  preprocessor.) 

Host  of  the  processor  cost  lies  in  basic  I/O;  by  integrating  the 
aacro  processor  with  the  language  translator,  this  cost  has  been 
ainiaized.  FKantorowitz  1976  ] Much  of  the  two-aan-aonths  spent 
in  iapleaentation  were  spent  in  understanding  nooks  and  crannies 
of  PL/I. 
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figure  2 


T is  not  intended  to  replace  any  existing  languages.  For 
distributing  mathematical  software,  Fortran  remains  the  only 
practical  medium;  for  character  processing,  something  like  PL/I 
or  SNOBOL  should  be  used.  Still,  for  the  bulk  of  scientific 
computation,  T ought  to  be  the  easiest  to  use,  particularly  since 
it  coexists  comfortably  with  Fortran  and  PI/I.  On  the  other 
hand,  cne  can  imagine  ways  that  T might  be  improved,  as  well. 

Features  omitted  for  ease  of  implementation  include: 

- trimmed  arrays,  like  X(2:N) 

- procedure  results  of  general  type 

- conditional  boolean  operators  that  dc  not  evaluate  their 

arguments  when  it  is  possible  to  avoid  doing  so 

- a swap  operator  * 

Fcr  other  features,  no  entirely  satisfying  design  was  apparent: 

- strings 

- more  general  procedure  calls  (such  as  indefinite  number  and 

type  o f arguments) 

- a means  of  constructing  arrays  directly  from  components,  as 
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a string  constant  constructs  a string  fros  individual  characters 
- a aeans  of  specifying  the  invocation  graph  of  who  calls  vhoa 

Perhaps  the  sost  fandanental  though  unavoidable  flaw  is  that, 
unlike  LISP,  the  language  is  not  trivial,  and  therefore  proqraas 
cannot  be  trivially  aanipulated. 


2.0  JUSTIFICATION  FOF  STILL  ANOTHER  SET  CF  GRAPHICS  PHIHITIVES. 
The  next  eiaaple  of  restyling  is  a staple  tot  reasonably  conplete 
interface  for  noninter  act  Ive  device-independent  graphics.  In 
addition  to  the  basic  line  drawing  primitives,  higher  level 
procedures  are  provided  for  displaying  functions  of  one  or  two 
variables.  This  interface  has  been  ispleaented  as  a library  of 
PL/I  procedures  which  call  the  SLAC  Unified  Graphics  package 
written  by  Robert  Beach  f 1978], 

Unified  Graphics,  with  its  eaphasis  on  the  ability  to  drive 
displays  like  the  IB  Pi  2250,  is  troublesoae  to  use  directly  for 
function  plots  and  the  like.  In  contrast.  Top  Drawer,  another 
graphics  systea  at  SLAC,  allows  for  function  plots  but  little 
else.  The  collection  described  in  detail  in  Appendix  G is  aeant 
to  strike  a useful  balance  between  these  twc  extremes,  and 
contains  aost  of  the  features  of  DISSPIA  important  for  scientific 
coaputation. 


2.1  ESTABLISHING  THE  ENT IRONHE NT . The  following  excerpt  froa 
the  example  prograa  given  in  section  1.1  akove  illustrates 
typical  preparation  for  plotting: 

STRUCTURE:  PF  * PLOT  FILE 

I NTEGE  R (500)  : HORN 

REAL  (2):  LL,  UR,  * CORNERS  OF  PECTA IGULAR  DOMAIN 

ORIGIN,  • FOCAL  POINT  FOR  SURFACE  PLOT 

X0,  SCALE  • COORDINATE  TRAN SFORH ATION  PARAHETERS 
G OPEN  (*  VE  P12FF'  ,PF) 

GPICT  (PF) 

SCALE  :-  0. 3333 
10  :«  -0.  5/SCALE  ( 1) 

3TRAN1  ( X 0 , SC  ALE , P F) 

The  plot  area  PF  is  used  to  reaeaker  various  options  and  to 
buffer  low  level  plotter  instructions.  This  work  area  is 
initialized  by  the  GOPEN  call,  which  specifies  the  output  device. 
(In  the  current  iapleaentation , no  corresponding  JCL  changes  are 
necessary.)  The  ease  with  which  devices  nay  be  changed  is  very 
useful  in  tuning  a plot  for  publication. 

For  coapatibility  with  numerical  procedures,  REAL  variables  are 
in  full  precision,  not  short.  At  the  start  of  each  new  picture, 
which  night  be  a screenful  on  a CRT  or  an  8.5  by  11”  page  on  an 
electrostatic  plotter,  GPICT  is  called. 
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All  plotting  is  done  relative  to  a user  coordinate  system,  which 
is  specified  by  calling  | 

GTRA  N ( F,  PP  ) 

where  F is  the  one  of  a procedure  which,  when  called  in  the  fora 
F(  X,  W,  PP  ) 
with 

SEAL  (M  | : X N <*  1 0 

BEAL  (2)  : W 

will  aap  the  point  X in  user  coordinates  into  a point  U in  the 
unit  square  [ 0,  1]x[0,  1 ].  normally  W ( 1)  is  thought  of  as 
horizontal  and  9(2)  as  vertical.  By  extending  PP,  the  user  can 
pass  parameters  to  P.  Por  convenience,  the  default 
transformation  naps 

W SCALE  * ( X - X0  ) 


2.2  DRAWING,  DIMENSIONING,  AND  FUNCTION  GFAFHING.  The  basic 
drawing  commands  are  G J HM  P , GDP  AW,  and  GTEXT  for  drawing  lines 
and  adding  text.  If  a nonlinear  coordinate  system  has  been 
soecified,  GCRAW  produces  a piecewise  linear  approximation  to  the 
implied  curve. 

A procedure  GGPAPH  is  provided  which  automatically  samples 
function  values,  sets  up  an  appropriate  scalinq,  graphs  the 
function,  and  dimensions  the  graph  using  round  numbers  in  a style 
consistent  with  the  format  used  by  T.  Pigure  3,  taken  from  Chan 
[ 1978  ],  is  a typical  plot. 

The  scheme  for  choosing  round  numbers  is  based  on  the  algorithm 
by  Dixon+Krormal  [ 1965  ],  Experience  and  an  informal  survey  of 
what  people  would  accept  as  being  "round  numbers"  led  to  various 
refinements.  As  in  Unified  Graphics,  the  choice  is  optimized 
over  a reasonable  number  of  major  tick  marks.  The  total  number 
of  tic  marks,  malor  and  minor,  is  not  allowed  to  be  either  too 
dense  or  too  sparse.  Por  a while,  the  number  of  minor  tick  marks 
was  chosen  so  that  each  interval  had  length  10**k,  but  for  input 
data  limits  (20,70)  the  resulting  tick  marks  were  at 
(-  130  ,0,  100,  200),  so  this  rule  had  to  be  relaxed  to  "either 
length  10**k  or  midpoint  of  major  interval."  If  the  difference 
between  the  data  limits  is  small  compared  to  the  magnitude  of  the 
limits  themselves  (as  occurs  for  example  in  plotting  a nearly 
constant  function),  then  the  labels  nay  tecone  unreasonably 
large.  Special  provision  is  made  for  th*s  case. 

Other  routines  are  available  for  scatter,  surface,  and  contour 
plots.  The  aontour  computation  uses  piecewise  quadratic  surface 
fitting  to  ensure  smooth  contours  and  proper  representation  of 
critical  points  ( "arlowvPowell  1976].  Pigure  « presents  output 
from  the  example  program,  which  computes  hill-climbing 
trajectories  for  a thr ee-gaussia n-peak  terrain. 
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CONCLUSION,  with  a level  of  effort  comparable  to  writing  a 
Fortran  preprocessor*  we  have  created*  by  conpiling  into  Pl/I*  a 
language  substantially  better  than  Fortran  or  its  derivatives. 
Since  PI/I  problees  cannot  be  altogether  avoided  by  this 
approach*  further  work  on  a language  like  T could  be  useful. 
Perhaps  the  effort  would  be  better  spent  cn  asking  LISP  a 
practical  language  for  scientific  coapuation  by  building  on  the 
research  in  syabolic  computation. 

Like  PL/I*  Unified  Graphics  is  good  for  a wide  range  of 
applications.  But  in  practice*  aany  people  won't  use  either. 

For  languages,  they  stick  to  Fortran;  for  graphics*  they  plot  by 
hani  or  not  at  all.  In  both  cases  it  has  proven  possible  to 
cheaply  restyle  the  existing  systen*  via  a preprocessing  phase  or 
driver  routines*  in  order  to  create  sore  agreeable  tools. 
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APPEND! I T 


Beport  on  the  programing  language  T 


"trie:  free  froa  anything  eatraneous; 

having  clean  lines  or  proper  proportion; 
the  state  of  readiness  for  action  or  use;" 

Webster's  Third  New  International 

"Everything  should  be  as  staple  as  possible, 
but  no  stapler." 

Einst  ein 

"What  it  lies  tn  oar  power  to  do, 

it  lies  in  our  power  not  to  do." 

Aristotle 

"In  all  spheres,  the  true  craftsaan  is 

the  one  who  thoroughly  understands  his  tools." 

Hoare 


TOKENS.  Prograa  text  is  nade  up  of  the  following  tokens: 


keyword  - one  of  the  following: 
CASE 
C01PLEX 
ELSE 
PIBST 
POP 

FORTRAN  PROCEDURE 

POBTBAN  PBOCEDUR E VARIABLE 

GET 

GET  ARRAY 
GET  DATA 
INTEGER 


NO 

PROCEIOPE 

FPOCEDURE  VARIABLE 
PUT 

POT  ARRAY 
POT  DATA 
FUT  DATA  ARRAY 
BEAL 

STROCTORE 

WHILE 

YES 


identifier  - a letter  optionally  followed  by  aore  letters 
and  digits. 

integer-constant  - one  or  aore  digits 

real-constant  - one  or  aore  digits  , a ".",  possibly  aore 
digits,  a "(",  pcssibly  a one  or  aore  digits,  and 

a ")".  Either  the  deciaal  point  and  succeeding  digits 
or  the  parentheses  and  exponent,  bat  not  both,  nay  be 
osltted.  Thus  1.,  0.23,  6.22<-23),  1 (-6)  , and 
3.  14159  (♦•001  wre  all  lsgnl  real-constants. 

string-constant  - a sequence  of  characters  enclosed  in 

apostrophes.  Apostrophes  are  not  allowed  in  the  string 
proper . 
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delimiter  - one  of  the  following: 

eoc  ()$*:.  , i | * *•/:■•  — <<■>>« 

(The  last  six  are  called  re  laticna  1-operators.)  4 

except  in  string-constants,  blanks  are  insignificant  and  nay  be 
used  freely  for  clarity.  Text  fret  a "•"  op  to  the  next  "eor" 

(end-of-record)  is  treated  as  a consent. 


PROCEDURES. 

prog  ran  * (procedure) 

procedure  ■ identifier  f " ("  identifiers  ")  " f result]  ] "eor* 
begin 

(declarat ion) 
f staten ent) 
end 

identifiers  * (identifier,)  identifier 
result  x identifier 
begin  « " ( ("  "eor"  . 

end  * ")  ) " "eor" 


Parameters  are  passed  using  call- fcy-r ef ercnce.  "result"  nay  be 
usel  just  like  any  other  for  Ml  paraneter,  in  particular  as  a 
destination,  but  nust  be  a scalar. 

?f  the  paraneter  list  is  oaitted,  the  procedure  is  assuaed  to  be 
the  top  level  aain  prograa. 


v 


DECLARATIONS. 

declaration  * scalar-type  [ • ("  bounds  ")"  ] ":"  identifiers  "eor" 
I "STRUCTURE"  ":»  identifiers  "eor" 
begin 

(declarat ion) 
end 


scalar-type 


I ["ECRTRAN")  "PROCEDURE"  [ " VARIABLE  " ) 

("()"  scalar-type  ] identifiers  "eor" 
* "INTEGER" 

| "REAL" 

( "CORPLEX" 


bounds  « [bounds,]  r expression : ] [expression]  . 


No  global  variables  are  allowed;  ooaaunicat ion  occurs  only 
through  paraneter  lists.  Ceclarations  reserve  storage  on  a 
stark;  the  variables  are  undefined  until  first  assigned  to. 


If  the  <ex pression>:  part  of  a bound  is  oaitted,  1:  is  assused. 
The  second  <expression>  should  be  oeitted  for  a foraal  array 
paraneter,  and  an  appropriate  value  will  be  taken  fron  the  length 
of  the  corresponding  actual  array  arguaent.  The  origin  of  an 
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actual  array  arguaent  need  not  natch  the  origin  of  th« 
corresponding  fornal  array  paraaeter.  For  esanple,  if  tha  actual 
arguaent  A was  declared  REAL  (0:7):  A and  the  fornal  paraaeter  B 
was  declared  B CAL  () : B,  then  B (i)  will  correspond  to  A (7). 

STRUCTU1ES  are  data  areas  assnaed  to  te  decoaposed  as  indicated 
ir.  the  subdeclarations.  Thus  if  actual  parameter  A is  declared 
STRUCTURE:  A 

IRTEGER:  I,  If 
REAL:  r 
couplet:  z 

and  corresponding  foraal  paraaeter  P is  declared 
STRUCTURE:  P 

IRTEGFE:  J,  L 
C3HPLEX:  H 

then  A. I and  P.J  correspond,  but  A.Z  and  P.V  do  not.  (If  F.  V is 
assigned  to,  both  A.X  and  A.Z  say  be  destroyed.) 

A procedure  VARIABLE  is  a variable  that  aay  refer  to  various 
actual  procedures;  in  contrast,  a PROCEDURE  is  literally  the 
naae  of  a procedure. 


STATE  RE  RTS . 

state  sent  ■ procedure-call  "eor" 

I destination  expression  "eor" 

I "GET("  identifiers  ")  " 

| "GET  ARRAY  <"  identifer  «)  " 

| "GET  DATA  ("  identifiers 
| " PUT  ( " arquaents  ")  " 

I "PUT  ARBAT  ("  destination  ") " 

I "PUT  DATA  ("  arguaents  ")  ■ 

| "RUT  DATA  ARRAY  ("  destination  ")" 

I "CASE"  "eor" 
begin 

(boolean-expression  "eor" 
begin 

(stateaent) 
end  ) 

("ELSE"  "eor" 
begin 

(stateaent } 
end  ] 
end 

| " WHILE  ( " boolean-expression  ")  * "eor" 

begin 

(states  ant) 

("FIRST"  "eor" 

(stateaent)  ) 
an  d 

( "FOB  (• 

((  expression  ("<"|"<«")  destination  ("<"("<■*) 
l(  expression  (">•(">*")  destination  (">"(">»") 


expression) 
expression)  ) "eor" 
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begin 

(statement) 

end 

GET  reads  from  the  next  record  of  the  input  data  a sequence  cf 
constants,  separated  by  cosaas.  For  GST  DATA,  each  value  should 
be  prefixed  with  "identifier  the  input  values  need  not 

appear  in  the  saae  order  as  the  corresponding  identifiers  in  the 
g FT  DATA,  and  if  a value  is  onltted,  the  variable  is  left 
unchanged.  POT  TATA  and  POT  write  out  the  current  values  of  the 
identifiers,  labelled  or  unlabelled,  in  an  intelligent  fashion. 

In  other  languages, 

CASE 
condl 
case  1 
con  12 
case  2 
ELSE 

easel 

light  be  vri  tten  as: 

IF  ( condl  ) THEM 
easel 

ELSE  IF  ( cond2  ) THE* 
ca  se2 
ELSE 

cas  e3 


Simla  rlv, 

while  cond 
part  a 
FIRST 
part  b 

night  be  translated  as: 

GOTO  FIFST 
TOP:  part  a 
FIRST:  part  b 

IF  ( cond  l THEE  GOTO  TOP 

If  the  FIFST  line  is  onltted.  It  is  assuaed  to  be  at  the  end  of 

the  loop;  that  is,  part  b is  eapty. 

Finally, 

EOF ( LOU  <-  I <•  HIGH  ) 
loop 

would  be  translated  as 

?OR  I * LOW  TO  HIGH 
loop 

and 

FOM  HIGH  > I >•  LON  ) 
loop 

as 

FOR  I • HIGH -1  BY  *1  TO  LOW 
loop 
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The  destination  and  expressions  controlling  a POP  loop  nay  not  be 
aodified  Inside  the  loop. 


expressions. 

destination  * (identifier  identifier 

) identifier  "("  subscripts  ") *• 

I i procedar  e-ident  if ier  . 
subscripts  * (subscripts,]  [expression] 

procedure- call  * procedure-identifier  ( " arguments  ") " ] . 

procedure-identifier  * identifier 

arguaents  * [arguments,  ] ( expression  | str  ing-const ant 

1 YES  | NO  ) . 

expression  * art  thaet  ic-  express  ion 
| boclean-expression 
a ri  thseti  c-expr  ess  ion  * f "-"  ] tern 
I a ri  thee  ti  c-  express  ion  ten 

I arithne tic-expression  tern 

tern  = factor 

I tern  factor 

I tarn  "/"  factor  . 
factor  = prinary 

I prinary  "•*"  prinary 
I prinary  prinary 

prinary  * integer-constant 
t real-constant 
| destination 
| procedure-call 
| " arithmetic-expression  ")  " 

boolean-expression  * r boolean-expr  ession  "|"]  boolean-factor  . 
boolean-factor  * (boolean-factor  "6"]  toolean-secondary  . 
boolean- secondar  y * ["-*"]  boolean-secondary 

I boolean-pri nar y 
bool e an-pr in  ar  y = YES  | NO 
| destination 
| procedure-call 

I ar  it  haetic-expression  relational-operator  ar  itheetic-expression  . 

If  a subscript  is  enpty,  it  is  assuaed  that  an  entire  row  or 
cclunn  is  being  referenced.  If  all  the  subscripts  are  enpty,  the 
parentheses  and  connas  nay  also  be  oaltted.  Array  expressions 
are  perforaed  eleaentwise. 

To  refer  to  a procedure  without  actually  invoking  it,  put  a * 
before  the  procedure  identifier. 

If  the  operands  are  nixed  INTEGER  and  REAL  the  result  is  REAL;  if 
either  is  COMPLEX,  the  result  is  COMPLEX.  Dividing  an  XRTBGCt  by 
an  INTEGER  and  raising  an  INTEGER  to  a power  are  illegal. 
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ppimitives  AMD  MACROS.  The  following  lacrcs  are  predefined: 
LENGTH  (arra  y,i)  size  of  array  in  the  i-tb  ditension 
DADD,  DHOL,  DDIV  eztended  precision  arithaetic 
CEIL,  FL03B  , SIGN,  ABS,  PLOAT,  RE,  IH 
“AX,  HIM  ( *1,  x2 , *3,  ...  ) 

POD  ( xl,  x 2 ) saallest  nonnegative  r such  that 
( «1  - r | / x2  is  integral 
PI,  EPS,  NAXRF.AL,  HAXINTEGEB 
AC3S,  ASIM,  ATAN,  COS,  EXP,  LW  or  LOG, 

LOG  10,  SIN,  SORT,  TAN 


The  following  procedures  are  predefined: 


NEXT  LIN E[  (n)  ] 


NEXT  PAGE 
END  CP  INPUT 

BLANK  SEPARATION  (n ) 

INTEGER  DIGITS  (n) 

PEAL  DIGITS  (n) 


PEAL  LEADING  DIGITS  f n) 


skip  to  start  of  next  line, 
then  put  out  n- 1 Jlank  lines 
(if  no  argusent  then  n*  1) 
pa  ge  e lect 

BOOLEAN  expression  set  to  TES  when 
the  input  file  tecoaes  eapty 
nuaber  of  blanks  to  be  left 
between  output  values  (default  3) 
nuaber  of  digits  for  integer 
output  (default  12)  1 <*  n <*  19 

nuaber  of  significant  digits  for 
real  output  (default  1)  3 <*  n <*  15 


0 <*  n (default  3) 


DATE  AND  TIME 


is  replaced  by:  'dd  aaa  19yy 


hh:  i 


Macros  are  evaluated  auch  like  procedure  calls,  first,  each 
arguaent  is  evaluated;  then  the  aacro  is  applied  to  ita  arguaents 
by  inline  expansion  according  to  the  aacro  definition;  finally, 
the  replaceaent  text  is  regarded  as  fresh  input.  During  the 
evaluation,  tokens  beginning  with  have  the  first  ■_ " stripped 

off.  (This  allows  aacros  to  be  teaporarily  "hidden.") 


DEPINE  ( identifier,  replaceaent  text  ) defines  a aacro.  The 
replaceaent  text  is  a sequence  of  tokens,  possibly  containing  eor 
and  aatched  parentheses.  Places  in  the  replaceaent  text  where 
arguaents  are  to  be  inserted  during  expansion  are  indicated  by  S 
followed  by  a digit  or  letter. 

IPELS  E (a,  b,  c,  d)  is  replaced  by  c if  the  token  a is  the  sane  as  b, 
ctharwise  by  d.  Either  a or  b aay  be  eapty.  for  exaaple,  the 
text 

DBPINE  (VERBOSE, TES) 

IPELSE  (VERBOSE,  TES,  PtlT  DATA  (X)  , ) 
would  be  replaced  by 
PUT  DA1A(X) 


A aacro  defined  before  the  first  procedure  applies  globally; 
other  aacros,  which  aay  teaporarily  redefine  the  global  ones, 
apply  only  to  the  procedure  in  which  they  are  defined. 
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APP  END  IT  I 


T implementation  notes 


FILE  PORHAT.  The  Report  Mentions  only  an  abstract  end-o f-record 
delimiter  and  not  any  particular  file  foraat.  This  is  because 
the  assuaptions  of  local  text  editors  are  of  overwhelming 
importance. 

The  current  implementation  provides  for  card  iaage  files,  in 
which  coluans  /3  through  80  are  i^iored.  In  order  to  obtain  text 
records  longer  than  72  characters,  continuation  lines,  flagged  by 
a blank  in  coluan  1,  nay  be  used.  Consents,  starting  with  a "t" 
end  in  coluan  7 2. 

Mora  technically,  each  card  has  an  end-of-line  character  SOL 
appended  to  it  and  each  record  has  an  BCF  appended  after  the  last 
FOL.  Thus 

REAL:  • constants 

PI,  EPS 

is  translated  to 

PEAL:  (^constants  EOL  PI,  EPS  EC  l EOF 

The  T processor  discards  blanks  and  text  fron  a • up  to  the  next 
ECL.  The  runtiae  I/O  package  discards  blacks,  coaaents,  and 
begins  and  ends.  On  the  other  hand,  the  PRINT  prograa  nakes  use 
cf  EOls  to  intelligently  fcraat  its  listing. 

’’’he  input/output  procedures  effectively  append  an  infinite  string 
of  EOF  characters  to  the  end  of  the  file.  The  characters  EOL, 

EOR , and  EOF  are  represented  internally  as  the  ASCII  control 
characters  US,  PS,  and  PS. 


OTHER  REHARKS  ON  USING  T.  Strings  eay  be  passed  through  a 
procedure,  for  exaaple  froa  a high  level  routine  to  a graphics 
trial  tive,  by  declaring  a foraal  paraater  as  INTEGER  (1). 

Various  options  aay  be  invoked  by  a sacro  call  of  the  fore: 
option  ( ? ) 

where 

? is  either  ON  or  OFF 

and  option  is  [default  value  is  in  brackets] 

RECURSIVE  [OFF]  recursive  procedures 

SUBCHECK  [OFF]  subscript  checking 

SHORT  [OFF]  short  precision  f 

FORTRAN  FACADE  [OFF]  procedure  looks  to  the 

outside  world  like  Fortran 

UNDERFLOW  [on]  turn  off  underflow  error  sesssges 

in  the  current  procedure  and  its  descendants 
(reals  that  underflow  are  set  to  0] 
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Arrays  say  have  at  aost  5 subscripts 


Por  the  purposes  of  STRUCTURE  all  eras  cat,  the  sizes  of  the  scalar 
1 types  are: 

INTEGER  1 unit 

REAL  2 units 

COMPLEX  4 units 

REALs  and  CCHPLEXe8  should  start  an  even  nuaber  of  units  into  the 

STRUCTURE  for  fastest  access. 


JCL  AT  SLAC.  In  order  to  invoke  the  PRINT  prograa  (designed  to 
intelligently  list  a set  of  T procedures),  use  the  PRINT 
catalogued  procedure  stored  in  WTL. CG.  EHG .E . An  ezaaple  of  a 
PRINT  run: 

//  JOB 

//PROCLIB  DD  DSN*NYL.CG .EHG.E, DISP*SHF 
//P  EXEC  PRINT 
//INPUT  DD  * 

...  T procedures  ... 

The  PRINT  prograa  recognizes  three  coaaands,  which  look  like 
cone  nts: 

•EE  causes  the  next  line  to  appear  on  the 

next  page 

•%T  title  causes  the  characters  laaediately  appearing 

after  "iET"  to  appear  in  the  title  line, 
and  the  next  line  will  appear  on  the  next 
page 

• XLn  causes  only  ideations  of  at  aost  level  a to 

appear,  e.g.  #1 (10  will  cause  indention  to 
be  suppressed 

In  order  to  invoke  the  precoapiler  and  coaplle  the  resulting 
PL/I,  use  the  TC  catalogued  procedure  stored  in  HYL . CG.  EHG. E.  An 
exaaple  of  a conpi le-and-go  run: 

//  JOB 

//PROCLIB  DD  DSN*NYL.CG .EHG.E, DISP*SHR 
//C  EXEC  TC 
//INPUT  DD  * 

...  T procedures  . . . 

//G  EXEC  PLIG 

Provision  has  been  aade  for  a PL/I  duap.  1o  get  this,  add  the 
JCL  card: 

//PLIDUNP  DD  SY SOUT*  A,  DCP*  (B ECPH*FB  A , LRECL*  13 3,BLKSIIE*  1 330) 
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APPENDIX  G 


Report  on  the  graphics  interface  G 


"In  improving  exploratory  data  analysis,  we  need  to  find 
new  questions  to  ask  of  the  data  (probably  the  hardest 
task),  and  new  ways  to  ask  old  questions,  throughout, 
arithaetic  as  a basis  for  preparinq  pictures  is  likely 
to  be  the  keynote.  It  is  sost  important  that  we  see  in 
the  data  those  things  we  do  not  expect  — pictures  help 
us  in  this  far  sore  than  numbers,  though  we  can  gain  a 
lot  lust  by  what  numbers  we  use." 

John  Tukey 


ESTABLISHING  THE  ENVIPONHENT.  In  order  to  remember  various 
options  and  settings,  such  as  character  sixe  and  line  type,  and 
as  a buffer  for  lo  wlevel  plotting  commands,  a work  area  PLOT  is 
provided  to  the  graphics  procedures.  This  say  be  declared  as: 
STRUCTURE  : PLOT 

INTEGER(SOO)  : NORN 
STRUCTURE:  USER 

e • e 

The  workspace  USER,  which  may  be  as  large  or  small  as  desired, 
allows  parameters  to  be  passed  to  user  procedures  called  by  G. 

To  initialize  PLOT  at  the  start  of  a run,  call 
GOPEN ( DEVICE,  PLOT  ) 
where 

DEVICE  is  a string  ccntainlng  one  of  the  codes: 

CALPICH  microfiche 

PDShO  1 3 Tektronix,  Hewlett -Packard 

VEP12PF  Versatec 

Note  that  if  VEP12FP,EXTS0RT  is  specified, 
then  the  reordering  of  the  plot  coamnnds 
necessary  for  the  Versatec  will  not  be  done 
by  a main  aenory  sort  (greatly  reducing  the 
run-time  aeaory  requirements).  In  this  case 
an  external  sort  step  must  be  supplied. 

Before  each  new  picture.  Including  the  first,  call 
GP ICT  ( Plot  ) 

Finally,  at  the  end  of  the  run,  call 
GC LOSE  ( PLOT  ) 

Omitting  this  aay  cause  the  last  picture  to  be  lost. 

Per  convenience,  the  initial  tra naf oraation  is  GTRAN1A,  which 
perforaa  the  sapping 

H :•  SCALE  * ( X - XO  ) 

whe  is 
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which  ace 


REAL  (2)  : H,  X,  XO,  SCALE 

and  XO  * (0,0),  SCALE  * (1,1)*  Tc  change  these  values, 
saved  in  PLOT.  WORK,  call 

GTBAN1(  XO,  SCALE,  PLOT  ) 

Line  types,  character  sizes,  and  character  angles  are  specified 
bv  calling 

GLTYPE(  LTYPE,  PLOT  ) 

GCSIZE ( CSIZE,  PLOT  ) 

GCANGL  ( CANGL,  PLOT  ) 
where 

LTYFB  is  a string  containing  one  of  the  codes: 

SOLID,  DOT,  DASH,  or  DOT-DASH 
REAL:  CSIZE,  • character  spicing;  initially  1; 

CAM  CL  • character  angle,  in  radians 

• counterclockwise  froa  horizontal; 

• initially  0; 


BASIC  drawing.  To  jusp  straight  to  the  pcint  X, 

G JUH P ( X,  PLOT  ) 
where 

REALM:  X • destination  (in  user  coordinates) 

To  draw  a line  fros  the  current  position  tc  X, 

GDFAWf  X,  FLOT  ) 

If  the  coordinate  transfor nation  is  curvilinear,  this  produces  a 
piecewise  linear  approximation  to  the  ivplied  curve.  Following  a 
change  of  coordinate  systes , GJUHF  should  be  called  before  GDRAW. 

To  write  out  text,  call 

GIEXI  ( X,  PRI,  SFC , PLOT  ) 
where 

REA L ()  : X • location  for  center  of  first  character 

PRI,  SEC  are  strings  of  length  at  nost  255 

In  order  to  obtain  a large  alphabet,  text  is  presented  to  GT EXT 
using  a pair  of  strings.  Every  pair  of  corresponding  characters 
in  the  priaary  and  secondary  strings  denotes  one  character  in  the 
extended  alphabet. 


The  secondary  character  for 

lower  Rob  an  upper  Ronan  lower  Greek  upper  Greek 
is 

L G H 

Por  cosnon  special  characters  the  secondary  character  is  also  a 
space.  The  priaary  character  for  Greek  letters  is  the  first 
letter  of  its  English  nase  or  one  of  the  special  cases: 

H eta  F phi  w oaega 

Q theta  T psi 

Additional  special  and  control  characters  are  awallable: 
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OC  begin  superscript 

1 end 

2 begin  subscript 

3 end 

4 save  position  1 

5 restore 

6 save  posit  ion2 

7 restore 

8 save  position3 

9 restore 

E increase  s ize 
P decrease 

IV  half  up 

2 down 

3 third  up 

4 dow  n 

5 sixth  up 

6 do  wn 

OU  one  Dack 

1 half  forward 

2 back 

3 third  forward 

4 bac  k 

5 sixth  forward 

6 back 

MX  is  an  cleaent  of 
H is  not  an  cleaent  or 
E there  exists 
\ for  all 
I intersect 
u union 

< is  strictly  contained  in 
> strictly  contains 
L is  contained  in 
B contains 

Uk  up  arrow 
0 down 
L left 
B right 

8 bidirectional 


IS  integral 
J contour  integral 
P partial-d 
D del 

5 plus-or-ainus 
X t ia  es 

: divided  by 

♦ abstract  ulus 

* abstract  tiaes 
2 radical 

0  infinity 
/ back  slash 
( left  sguare  bracket 
) right  square  bracket 
L left  anqle  bracket 
R right  angle  bracket 
k left  curly  bracket 
Z right  curly  bracket 
< less  or  equal 
■ not  equal 
E equivalent  to 
Q proportional  to 
> greater  equal 
0 degree 
E section 

6 dagger 

P double  dagger 
H h bar 
V last  da  bar 

0 underscore 
v overscore 

OP  cross 

1 diagonal  cross 

2 diaaoed 

3 box 

4 star 

5 diaqonal  start 

6 cross  with  serifs 

7 diaqonal  cross  with  serifs 

8 coapass  rose 

9 octagon 


« 


For  exaaplc,  the  definitioc  of  the  gaaaa  function  is  given  by: 

GTEXT(X,  • (R-1)!*G  (E ) -T  400  52F  0*  E O-TITOH-1 IDT* , 

• L B L SCCSCCC  C IC  ICICI  CLL*,PL0T)  * 
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DIMENSIONING.  Given  Halts  on  the  data  canoe,  suitable  values 
fcr  plotting  the  data  saf  be  obtained  by  calling 

G SCALE  ( DATA  HIM,  DATA  HA  X,  # input 

LABEL  HIM,  LABEL  HAX,  EXP,  HAJCB,  HI  HOB  ) • output 

The  data  Mill  then  be  bracketed  by  LABEL  HIR*K)**EXP  and  LABEL 
HAX  * 1 0**EXP,  which  ace  round  nusbers  in  engineering  fornat. 

GTU  acts  3oaewhat  like  a ruler,  drawing  an  unlabelled  axis  with 
large  and  snail  tic  aarks: 

GTIC  ( L,  H,  HAJCB,  HIHOR,  OFF,  PLOT  ) 
where 

BEALO  : L,  H,  • endpoints  of  axis; 

CFF  • offset  coordinates  of  najor  tic  sack 

• endpoints,  which  deternine  the  sixe 

• and  direction  of  the  tic  aarks; 

BEAL:  NAJOR,  HIMOR  • nuaber  of  aalor  and  ainor  tic  aarks, 
t counting  endpoints;  thus  a yardstick 

• aight  have  HAJ0B*4,  HTN0P*13; 

GT.A3  adds  integer  labels  to  the  tic  aarks  produced  by  GTIC: 

GLAB  ( L,  H,  NAJOR,  OFF,  LOU,  HIGH,  PICT  ) 
where 

INTEGER:  LOW,  HIGH  f label  values  at  endpoints; 

REALO  : OFF  • offset  of  first  character  of  label 

L,  H,  HA JOB  are  as  in  GTIC. 

GF0BH1  lays  out  a fora  suitable  for  scatter  plots  and  function 
qn  ph  s: 

GF0R.1 1 ( GPP  I,  GSEC, 

XPRI , XSEC,  XLOW,  XBIGH, 

TPPI , TSEC,  FLOW,  THIGH,  PLOT  ) 

where 

GPRI,  GSEC,  XPRI,  XSEC,  TPPI,  TSEC  are  pairs  of  strings 
defining  the  general  title,  X-  and  T-axis  labels; 
REAL:  XLOW, XHIGH, FLOW, THIGH  • specifies  the  data  liaits 
GFoanl  autoaatically  sets  up  a coordinate  systea  so  that  the  plot 
will  fill  the  screen.  To  add  a function  curve,  just  use  GJOHP 
and  GDRAW. 


PLOTTING.  G SCAT  provides  scatter  plots: 

GSCA T ( GPRI,  GSEC, 

XPRI,  XSEC,  X, 

TPR I , TSEC,  T,  PLOT  ) 

where 

REALO  : X,  T • data  points 

GPRI,  ...  , TSEC  are  as  in  GF0BH1. 

ggbaph  provides  graphs  of  functions  of  one  variable: 

GGBAPH  r GFPI  , GSEC, 

XPPI,  XSEC,  A,  B, 

FPRT , FSEC,  F,  PLOT  ) 

where 
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: 


BEAL:  A,  B • endpoints  of  interval  for  graphing 

FROCEDURE  ()  BEAL:  F • function  to  be  plotted,  which 

• is  called  in  the  fora: 
f T * F ( X,  PLOT  ) 


To  plot  contours  of  the  surface  passing  through  the  points 
( X (I)  # T <J)  , P (I,J)  ) 

c a 11 

GCON  T ( LL,  U8  , F,  LEVELS,  FLOT  ) 
where 

BEAL  (2):  LL,  UP  • coordinates  of  lower  loft  (aost 

• negative)  and  upper  right  (sost 

• positive)  corners  of  the  rectangle 

• in  the  X T plane  on  which  data 

• is  given 

B EAL  ( , ) : F • data  values 

REAL  ()  : LEVELS  • contour  levels  tc  be  plotted 
A uni  fora  grid  is  assuaed. 


Tc  draw  a transect  surface  plot  with  hidden  lines  reaoved,  call 
GSUB  F ( LL,  UP,  F, 

AZIMUTH,  ELEVATION , OPIGII,  SCALE,  PLOT  ) 

where 

BEAL:  AZIMUTH,  ELEVATION,  SCALE 
BEAL  (2)  : CP  I GIN 

The  coordinate  tra asf oraa tion  used  in  GSUBF  naps  P * (X,Y,P)  into 
ORIGIN  ♦ ( -SI  Cl  0 ) (P-C) ‘SCALE 

( -S2*C  1 -S  2*S  1 C2  ) 

where  Cl  = COS(AZIMUTH)  , ...  , S2  * SI N (ELEVATION) 

C * ( (LL  ( 1)  »UP  ( 1)  ) /2  , (LL  (2)  ♦UP  (2)  1/2,  0 ) 

LL,  UB , and  F are  as  in  GCONT. 
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